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Abstract 

The Klein-Gordon system describing three scalar particles without 
£^ ■ interaction is cast into a new form, by transformation of the momenta. 

Two redundant degrees of freedom are eliminated; we are left with a 
covariant equation for a reduced wave function with three-dimensional 
arguments. This new formulation of the mass-shell constraints is equiv- 
alent to the original KG system in a sector characterized by positivity 
QJ ] of the energies and, if the mass differences are not too large, by a mod- 

erately relativistic regime. 

Introducing mutual interactions provides a model which is (at least for 
^ . three equal masses) tractable and admits a reasonable nonrelativistic 

H ' limit. 
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1 Introduction 

1.1 Motivations 

Relativistic particle dynamics is concerned by situations where the par- 
ticles we consider are not significantly created or anihilated, whereas 
other relativistic effects must be taken into account. In principle the 
description of such particles should result from a specialization of quan- 
tum field theory (QFT) to its n-body sector. 

This line leads to the famous integral equation of Bethe and Salpeter 
(BS) in the two-body case. Three-body generalizations have soon been 
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considered in the litterature [1]. More recently see [2] [3]. For n > 2 how- 
ever, the complexity of the BS approach seems to be almost prohibitive 
as far as practical applications are concerned. 

An alternative approach, based upon first principles [4] [5], uses n mass- 
shell constraints in the form of coupled wave equations where interaction 
terms can be either phenomenological or derived from QFT [6]. This 
method shares with BS equation the property of manifest relativistic 
invariance, realized at the price of dealing with redundant degrees of 
freedom, since the arguments of the wave function are four-vectors. In 
the two-body case, there is a clue for eliminating the redundant degree 
of freedom: the sum of wave equations rules the dynamics, whereas 
their difference allows to determine how the wave function depends on 
the " relative time" . This dependence turns out to be trivial and one is 
left with a three-dimensional problem. 

In the three-body case we have to cope with two "relative times" . These 
superfluous degrees of freedom are present as well in the three-body ver- 
sions of the BS equation. Their elimination (or factorization) is desirable 
for physical interpretation; it would produce (after diagonalization of the 
total linear momentum) a reduced wave equation which is covariant but 
similar to a Schroedinger equation with three-dimensional arguments. 
Unfortunately, the simple procedure utilized in the two-body case does 
not work for n > 2. 

An important issue of n-body dynamics is cluster separability; but a 
less restrictive and more essential requirement is global separability: 
one must at least recover free-particle motion when all interactions are 
put equal to zero. Models violating global separability have been con- 
sidered in the past [7] [8] , mainly for their computational simplicity, but 
we belive that any reasonable formulation of n-body dynamics must in- 
clude free motion as a limit when all the terms carrying interactions are 
"switched off'. 

In sofar as fermions are concerned, these matters have been discussed 
earlier in the literature [9] [10]. 

For scalar particles with masses m a , free motion can be described by n 
Klein-Gordon (KG) equations, say (p„ — m„) 3? = where <3? depends 
on the momenta p±, ....p n . We can give a sharp timelike value k a to the 
total linear momentum and use the differences of these equations. In 
the two-body case, it follows that the relative time (or alternatively the 
relative energy ^(pi — P2) ■ k/V~k?, which is conjugate to it) arises only 
in a trivial factor of the wave function. 

But this procedure is unable to produce any simplification as soon as n > 
2. So we face this difficulty that even for free particles, the usual form 
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of the equations of motion fails to permit the elimination of superfluous 
degrees of freedom. 

This point may seem to be academic, because a system of noninter- 
acting particles has no bound state, which renders a three-dimensional 
formulation unnecessary. But we bear in mind the eventuality of intro- 
ducing interactions that ultimately give rise to bound states. Therefore 
the possibility of a reduction is essential and should survive in the free 
case. 

In this paper we focus on three-body systems and we firstly consider 
the case of noninteracting particles. Let us stress that the free system is 
not considered on its own right, but rather as preliminary to the further 
introduction of mutual interactions. 

Since the KG equations as they stand do not permit a factorization 
of the dependence on relative times, it is natural to transform these 
equations into an equivalent system such that two superfluous degrees 
of freedom can be desentangled from the kinematics. 
An early attempt to carry out this task for an arbitrary number of par- 
ticles was made by Sazdjian [10] [11] fifteen years ago. Here, however, we 
shall be concerned with the 3-body case only, and shall take advantage 
of a simplification that is not possible for n > 3. 

Our aim is to eliminate two degrees of freedom in such a way that 
the mass-shell constraints reduce to a covariant problem with three- 
dimensional arguments. Ultimate introduction of interactions will be 
briefly sketched at the end. Of course, the Poincare invariance of kine- 
matics must be preserved and all particles should be treated on equal 
footing (democracy). These conditions are not likely to select a unique 
scheme, but if we intend to make it as simple as possible, there are not 
too many choices. 

We perform a rearrangement of the individual coordinates (well known 
in celestial mechanics) which is adapted to the consideration of relative 
variables. We insist on having invertible formulas, which is necessary 
in order to make sure that the new form of the equations of motion is 
equivalent to the original KG system. 

Section 2 is devoted to an exposition of the notation used and of the 
basic useful equations of relativistic dynamics. In Section 3 we collect 
known results and perform elementary manipulations. 
In Section 4, using the "heliocentric variables", we construct in closed 
form a transformation of the free-particle system and discuss under 
which conditions this transformation is invertible. In Section 5, we 
briefly indicate how mutual interaction could be introduced. 
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2 Basic equations, notation 



Units are such that fi = 1 whereas c remains unspecified. 
We start from the KG-system describing n particles in momentum rep- 
resentation 

p 2 a § = m 2 c 2 § a,b,c=l...n (1) 

where depends on the three four- vectors p%. Configuration and mo- 
mentum variables are mutually conjugate [q%,Pb@] = ifidbSp, and so on. 
We make use of the following notation: 

Q=-^}Zla, P = J2p a , z ab = q a -q b (2) 

v. = (3) 

n 

Moreover it is convenient to define 

Pab =Pa+Pb Vab = ^ (Po ~ Pb) (4) 

Beware that z a b is not conjugate to y ab . We obviously have the following 
relations 

Va-Vb = ~{Pa - Pb) = -tyab (5) 
X] Va = °> 2 Vab + Va-yb = 0, Pa = ~Va 

The tilde symbol denotes projection orthogonal to P a , in other words 
y a = II y a , z a =Uz a , with II = 5 - (P ® P)/ P 2 . Similarly the "hat" 
symbol refers to the projection orthogonal to k a , eigenvalue of P a . For 
instance 

VI = y a a -(ya-P/P 2 ) P a 
r a =y a a - (y a • k/k 2 ) k° 

Heliocentric variables 

The problem of " relative times" cannot be easily handled unless we first 
choose a set of independent relative variables. For this end, one particle 
is arbitrarily picked up; let it be particle with label one. With respect 
to particle 1, the relative configuration variables are defined like in [7]. 

za = qi - qA (6) 

where the capital labels A, B, C run only from 2 to n. From (3) it follows 
that za is conjugate to yA- 
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Let us now specialize to three-body systems; we can write 

1 1 

2/12 = 2/2 + 7^2/3 2/13 = 2/3 + 7^2/2 (7) 

^12 + ^23 + ^31 = (8) 

Notice that eqs (1-5) hold true for any n, whereas (7) (8) are valid for 
n = 3 only. It is clear that Q, z 2 , 23 are independent configuration vari- 
ables. In the same way P, 2/2>2/3 are independent momentum variables, 
canonically conjugate to them. We can use the set of canonical variables 
Q,z 2 ,z 3 ,P,y 2 ,y3 in place of qi,q 2 ,q3,Pi,P2,P3, this change is trivial. 
In this "heliocentric" formulation, democracy among the three particles 
is of course not kept manifest but can be checked at various stages of 
the development. A similar re-arrangement, showing up two relative 
momenta, is of current use in (Newtonian) celestial mechanics. 

Among the quantities P a b we shall more specially need to evaluate 
P\2 , P\z . They are given by 

Pyi = \p + V3 Py3 = \p + V2 (9) 

We shall also need the canonical expression of 2/12,2/13) given by (7). 
It will be convenient to replace eqs (1) by their sum and their differences; 
to this end we define , 

V A = -(ml- m 2 A ) 

so the equal-mass case is characterized by the vanishing of both v±, v 2 . 



3 Equations of motion. 

Equations (1) can obviously be written 

C 2 E™« $ = Y.Pa * (10) 

(ml-ml)c 2 <$> = { P l-p 2 b )<5> (11) 
Notice that, according to notation (4) 

\{pl-pl)=yab-Pa b (12) 

In equation (10), let us use the identity 

n 

nY.P 2 ^P 2 + Y.(Pa- Pb ) 2 (13) 

1 a<b 
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valid for any sum of n squares in a commutative algebra. We obtain 

3^Tm 2 c 2 $ = P 2 ^ + Y.(Va-Vb? * (14) 

a<b 

In terms of the relative variables (see eq. (7)) we have another identity 
specific of the three-body problem 

J2(p*~Pb) 2 = Q(yi + vi + V2 ■ y 3 ) (15) 

a<b 

Now in the r.h.s. of (15) we separate time from space according to the 
direction of P, and insert the result into (14). We get 

£(p«, - Pb? = D + 6P 2 E (16) 

a<b 

where 

D = 6(y 2 + y\ + y 2 ■ y 3 ) (17) 

S = {P 2 r\y2 ■ P? + (ys ■ P? + 0/2 • P)(y 3 ■ P)\ (18) 
Thus the sum of eqs (1) is 

(3^m 2 c 2 -P 2 ) $ = (D + 6P 2 E) $ (19) 

The remaining combinations of (1) can easily be written as the "differ- 
ence equations" 

2/12 ■ P12 $ = ^2C 2 $, J/13 ■ i> 13 $ = z/ 3 c 2 $ (20) 

Now it is natural to require that $ is also eigenstate of the total mo- 
mentum, say 

P a q> = k a $ (21) 

for some timelike constant vector k. But (in contrast to what happens in 
the two-body case) this procedure is unable of getting rid of the relative 
energies cyA • k/Vk? conjugate to the relative times c" 1 za ■ k/Vk 2 ~. 

Nevertheless, we can look for a new set of canonical variables; if these 
variables are suitably choosen, equations (20) may after all result in the 
elimination of two degrees of freedom. 
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4 Alternative Formulation of the Free 
Motion 



4.1 Transformations in momentum space 

We shall construct a new representation of the KG system. It will in- 
volve a new set of operators q' a ,p' b satisfying the canonical commutation 
relations. Let them be rearranged as P',z' A ,y' B by formulas similar to 
(2) (3) (6). In particular v' a van ishes and P' = J2p' a but we must 
require that P' = P in order to preserve translation invariance. Thus 

P 1 

y'a = --Pa V'ab = ^{Pa-P'b) 

y'a-v'b = -{P'a-Pb) = - 2 y'ab 
3 

Naturally X) y' a = 0. Notice for V^p' and for Y^(p„ — Pb) 2 identities 

1 a<b 

similar to (13) and (15). We obtain 

y'u = V2 + \y'z via = y 3 + \y'i (22) 

4 2 4 2 

y'2 = 3Z/12 - 3^13 ^3 = 3-2/13 - 32/12 (23) 

Define 

Q' = ^E^ A = 4l-4a (24) 

It is clear that Q' , z' 2 , z' 3 , P,y 2 ,y 3 must be independent variables, 
conjugate to z' A , etc. 

A transformation in momentum space will be enough to induce the suit- 
able transformation among operators. In fact we are going to construct 
the quantum analog of a point transformation in momentum space (see 
Appendix 1). 

Let us start with a wave function &(pi,P2,P3)- Perform a change in the 
space of its arguments 

Pa ^ Pb 

or equivalently P 1— > P' = P and yA l— ► y' A - 

d d 
Instead of the old configuration variables z = i— , Q = i- 



we shall now consider 



'dy ^ dP 



z' = i— Q' = i— 

z l d y > * 8P> 
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Since dP/dy' = and dP/dP' = 8, the transformation formulas are as 
follows 

z'1 = zb. (25) 

Q' a = Q a + ^ z AiX (26) 

with summation also over (repeated) capital indices. In these formulas 
it is clear that the transformation of momenta must be invertible. 
Beware that Q' may not coincide with Q because of dy/dP'. In addition, 
we observe that the new relative coordinates actually mix the old ones. 
However, we shall prove later (Section 5) that this difficulty disappears 
in the large-total-mass limit. 

It is in order to stress that finding the desired transformation amounts 
to solve a problem in the framework of c-numbers. The question of 
inverting formulas, discussed below, is nothing but a nonlinear problem 
concerning the arguments of the wave function. Since it is specified that 
we are dealing with momentum representation, we shall use without 
confusion the same symbols for the arguments of the wave function and 
the multiplicative operators they define. 

For a better understanding of the mathematical structure, it is perhaps 
relevant to notice that Q' and z' A are "formally hermitian" in this sense 
that they are symmetric operators in 

£' 2 (R 12 ) = L 2 (R 12 , d*Pd% d%) 

whereas Q and zb are symmetric operators in 

£ 2 (R 12 ) = L 2 (R 12 , d 4 P A2 As) 

In contrast the momenta are symmetric operators in both senses. For 
mathematical convenience we shall work with a new wave function 

* = |J|2<3>, where J is the Jacobian J = ^TT^T^ll always finite 

D{Pi,p2,P 3 ) 

and nonvanishing insofar as our transformation is invertible. Indeed 
multiplication by \J\i maps C 2 onto £' 2 , so $ (resp. \£ ) belongs to 
the rigged-Hilbert space constructed by taking C 2 (resp. C' 2 ) as Hilbert 
space. Although C 2 (resp. C' 2 ) has no direct physical meaning, it allows 
for representing the Poincare algebra and gives a rigorous status to the 
operators involved in the wave equations. 

Since p a are multiplicative operators, they commute with J, so the mass- 
shell constraints can be written either as (1) for $ or equivalently in the 
form (p 2 — m 2 c 2 )^ = 0, with each p a expressed in terms of p' b . 
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In momentum space, the Lorentz group is characterized by this property 
that it leaves all the products p a ■ pb unchnged. 

Provided all the p' a • p' b can be expressed as functions of p c ■ pd and vice- 
versa, the same realization of the Lorentz group can be as well charac- 
terized by invariance of all the scalar products p c -Pd- in such a situation, 
although M' = J2l'^p' ma Y be distinct from M = Q A p, their com- 
ponents span the same Lie algebra. Moreover J being conserved by 
rotations, it follows that M and M' are both symmetric in C 2 and also 
in £ . 

Till now we have considered a large class of transformations, charac- 
terized by equations (25) (26); the classical (non-quantum) limit of such 
formulas would define point transformations in momentum space. 
We now specialize to a transformation which allows for eliminating the 
superfluous degrees of freedom. All we need is an invertible transforma- 
tion such that 

P12 ■ 2/12 = P ■ y[ 2 Pi3 -vn = P- y'i 3 (27) 
Indeed, if these relations are satisfied, (20) takes on the form 

y' 1A • P * = v A c 2 (28) 
Then according to (23) the "difference equations" are 

l/ 2 -P* = (^-|^)c 2 * (29) 

^.p* = (^ 3 -^ 2 )c 2 ^ (30) 
With help of equation (21) we obtain 

* = 

5(P a - k a ) 5(y' 2 ■ k - ^v 2 c 2 - 2 -vzc 2 ) S(y' 3 ■ k - K z c 2 - ^ 2 c 2 ) V (31) 

where ip depends on y^y's only through their orthogonal projections 
onto the three-plane orthogonal to k. One remains with the problem 
of determining a reduced (or internal) wave function ip which has no 
more arguments than the wave function of a nonrelativistic problem. 
The dependence of the wave function on y' A ■ k is now factorized out. 

For simplicity we complete our transformation law by imposing that 
the space projections of 2/2,1/3 (with respect to the rest frame) remain 
unchanged, say 

y' A = yA (32) 
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and just transform their time components in the way dictated by eqs 
(27). This choice obviously preserves Lorentz invariance; we shall prove 
below that it does not destroy the democracy among particles. 

In view of equation (32), and taking into account the identity 

VAa^^^Pa+VAa (33) 

it is clear that our change of variables is essentially determined by (27). 

As they stand, these formulas implicitly define y' 2 • P and y' 3 • P in terms 

of the old variables; but we still have to solve (27) for y 2 ■ P and y' 3 ■ P 

in order to exhibit the transformation in closed form. 

According to (9) and (7) the left-hand sides of conditions (27) are as 

follows: 

2 11 

P\2 ■ 2/12 = g2/2 • P + 2/2 ■ 2/3 + • P + (34) 

2 11 

^13 • 2/13 = ^2/3 • P + 2/2 ■ 2/3 + 7^2/2 • P + -jyl (35) 

For the right-hand sides, eqs (22) yield 

y'i2-P={y'2 + \y'z)-P (36) 
y'iz-P = i.y'z + \y'2)-P (37) 

Therefore the requirement that (27) are satisfied can be expressed as 
the linear system 

\v2 ■ p + V2 ■ ys + • p + \y\ = {y'2 + \y'3) ■ p (38) 
\y3 ■ P + V2 ■ V3 + \y2 ■ P + \y\ = (y 3 + \y'2) ■ P (39) 

to be solved for y' 2 • P and y' 3 • P. The outcome of system (38) (39) is 

y'2-P = \{y2-P + y2-V3 + yl)-\y 2 2 (40) 
y'3-P = \(yz-P + y2-y3 + yl)-\yi (41) 

whereto we insert the decomposition (33). This substitution, together 
with (32), determines in closed form the transformation of momenta. 
But it remains to be checked that this transformation is invertible. 
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Translation invariance was ensured from the outset by assuming that 
P' = P. 

Lorentz invariance is preserved because all the quadratic scalar quanti- 
ties formed with the vectors P, 1/2,2/3 are scalar invariant under space- 
time rotations. 

Democracy between particles is not manifest in the heliocentric nota- 
tion. Nevertheless it is not difficult to check that our way of transform- 
ing momentum variables treats all three particles on the same footing. 
Indeed we first observe that (32) entails y[ = y±, which amounts to 
finally write y' a = y a for the three particles. Then using (12) and (4) we 
realize that (27) automatically imply a third relation P23 ■ 2/23 = P • J/23- 



4.2 Inversion of formulas 

Now that all components of the new momenta p' a are determined we 
can (in principle) evaluate the configuration variables through formulas 
(25) (26). It is essential to realize that our transformation of the mo- 
menta among themselves must be invertible: if it were not, the trans- 
formation would not be canonical and the new form given to the wave 
equations would not be equivalent with the KG system. 
Formula (25) can be written in closed form provided we are able to carry 
out this inversion. We are thus faced with the problem of mapping the 
new momenta back onto the old ones, which amounts to solve the system 
(38) (39) now for the unknown 2/2 ■ P, 2/3 ■ P in terms of y' 2 • P, y' 3 ■ P, 
assuming this time that the latter are given and taking (32) into account. 

Positive-energy condition 

The domain where (38) (39) must be inverted can be limited to the 
positive-energy sector. So we require not only that P is timelike and 
future oriented, but also that every vector p a is timelike and points 
toward the future, which entails P ■ p a > and p a ■ Pb > 0. 

At this stage it is convenient to introduce the dimensionless quantities 
thus (18) becomes 

~=e + 7 ] 2 + Vt (43) 

The positive-energy condition above implies limitations for £ and rj. 
Indeed we first derive from (3) 

£ = I-^ V = ±-^* (44) 

? 3 P 2 1 3 P 2 { ' 
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From positivity of P • pa we get 



£ < I, V<1 (45) 



On the other hand we have P ■ pi = P 2 — P • P2 — P ■ P3 • According to 
(44) this identity reads P ■ pi/P 2 = — h£ + ?7 and this expression also 
must be positive. We end up with 

~1<Z + V (46) 

With these limitations in mind, we now turn to the inversion of system 
(38) (39). In view of the identities 

(y A • P)(y B -P) _ 
VA-VB = ^2 ^VA-VB 



we can write 



2/2 ■ 2/3 = P 2 £v + 2/2-2/3 



vl = P 2 e + vl vi = PW + vl 

Insert these formulas into (34) (35), and write (27). We get 

Because of (32) all quantities of the form yA ■ Db are already known. 
The above system (47) (48) is quadratic in the unknown quantities £,77. 
Define dimensionless quantities u, v through the formulas 

P 2 u = y / 12 -P-(y 2 -y 3 +^y 2 ), P 2 v = y' 13 -P- (2/2 ■ 2/3+^2/1) (49) 

They are regarded as functions of the new momenta, since P and yA 
coincide with P' and y' A respectively. Inserting (22) into (49) yields 

P 2 u = y' 2 ■ P + \y' z ■ P - y' 2 ■ y' 3 - ±(y> 3 ) 2 (50) 

P 2 v = y' 3 -P + \y' 2 -P-y' 2 -y',-l(y' 2 ) 2 (51) 

The system (47) (48) becomes 

|C + \v + £77 + *L = u (52) 
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2 1 £ 2 

3»/+3^ + ^ + y=« (53) 
to be solved for £, rj with u, u as in (49). Setting 

3(u + v) = a, S(u — v) = e 

system (52) (53) can be cast into the form 

( + V + 2<l( + ^(( 2 + V-) = | (54) 

St-^ + itf -<*) = ! («) 

It is convenient to define X = £ + 77, Y = £ — 77. When inserted into 
(18) this change of variables produces 

■~ = -X 2 + \y 2 (56) 
4 4 

System (54) (55) becomes 

\X 2 + X- Y = I (57) 

2Y - 3XY = 2e (58) 

The positive-energy conditions (45) (46) demand that X belongs to the 
1 2 

open interval (--, -) and also that Y > —1, which in turn require that 
*<§(l + e). 

2 

When e = a couple of obvious solutions is given by X = — (whatever is 

O 

a) which corresponds to Y = ±2yl — cr/3, but this possibility is ruled 
out by (45). Other solutions are given by Y = hence 

X = X ± = -(-1± v / TTct) (59) 

but the solution X - is excluded in view of condition (46). 

We now turn to the general case. The possibility that strictly X = 2/3 
being discarded, we now solve (58) 

Y = — ^ (60) 
2-3X v ; 
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and bring the result into (57). Hence a 4th-degree polynomial equation 
to solve for X, 

(2-3X) 2 (^X 2 + X-^) = e 2 (61) 

Graphic analysis 

In principle such equation can be explicitly solved by radicals. But a 
graphic analysis gives a better understanding. Solving (61) amounts to 
discuss how, in the X, Z plane, the parametrized curves Z = R n (X) = 

(2 — 3X) 2 (-X 2 + X — — ) are intersected by a straight line with trivial 

equation Z = e 2 . See Appendix 2. 
The outcome of graphic analysis is: 

Proposition I. 

3 1 11 

Provided — - < a < - and e is taken in the open interval (—-,-), 

among the real solutions of the system (57) (58) there exists a unique 

1 2 

one, X, Y such that X £ ( — , -) and such that X reduces to X + when 

3 3 

e vanishes. 

Moreover we observe that 

2 2 

X e remains bounded by —, ensuring that Y > —1 as required among 

3 3 
the positive-energy conditions. 

The expression X = S(a, e) for this solution could be written in closed 
form, but is very complicated, except naturally for vanishing e where it is 
just given by X + . For applications, we have better to use a development 
in powers of e 2 , say 

X = S(a, e)=X + + e 2 X (1) + e 4 X {2) + e 2 ^X {v) + .... (62) 

All coefficients X^ are derived from (61) and depend on a. We find 

4 

for instance X (1) = + _ _ g 2 . Note that 

S = - + 0(a 2 ,e 2 ,ea) (63) 

For the sake of a physical interpretation, investigating the behavior 
of our formulas at large P 2 is of interest. Equations (50) (51) show that, 
considered as functions of the independent variables y^ Q ,P^, all the 
quantities u, v, a, e, X, Y are of the order of 1/|P|. We simply have 

C = 2u-v + 0{l/P 2 ), n = 2v -u + 0{l/P 2 ) (64) 

Proposition (I) stated above ensures that the transformation from the 
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old momenta to the new ones is safely invertible in an open set of values 
given to the couple a, e. As these quantities are first integrals for free 
particles, their limitation to an interval defines a sector which is invari- 
ant by the motion. Characterization of this sector in terms of physical 
quantities will be discussed in the next Section. 
Remark: 

Infinitely many other domains ensuring a unique solution to (57) (58) 
could be exhibited. But we can enlarge the interval for e only at the 
price of shrinking the one for a. 

4.3 Physical conditions 

In view of eqs. (28), the wave function includes a factor 5(y' 12 ■ P — 
c 2 u 2 )8{y'i 3 • P — c 2 f3). The relevant domain for the arguments of \I/ is 
thus limited by the constraints y' 1A • P = (?va-, where the masses are 
given from the outset. On the mass shell, we can replace y' 1A P by c 2 ua 
in (49) or in the definitions of a, e. 

The particular case where e vanishes is interesting because it arises when 
the particles are mutually at rest, provided m 2 = 777.3, which includes 
the special case where all masses are equal. Moreover e remains small 
insofar as U2, ^3 and the velocities are not too large. 
For simplicity, let us focus on the assumption that u 2 , ^3 are small 
enough. In order to keep some contact with nonrelativistic mechanics, 
our scheme must encompass the case e = 0; thus the solution which 
reduces to X~ for vanishing e is excluded. Since the transformation of 
momenta must be one-to-one, we are also obliged to discard the solutions 
which reduce to the fixed point for vanishing e. 

Finally we have no other choice than the solution given by X = S(a, e). 

Let us now discuss in more details how we can manage, by simple 
physical requirements, to keep a, e within admissible values allowing to 
apply Proposition I. 
From (49) we obtain 

P 2 | = (u 2 + u 3 )c 2 - (2y 2 ■ y 3 + X -y\ + l -y\) (65) 

P 2 { = {v2-vi)c 2 + \{t 2 -yl) (66) 

in other words 

3 11 

a = a o ~ pi ( 2 ^ 2 " ^3 + + 2^3) ( 67 ) 
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e = £ o + ^ (yl-y!) (es) 

setting 

ao = 3(iA2 + ^ 3 )c 2 /^ 2 , eo = 3(^ 2 - ^ 3 )c 2 /^ 2 (69) 

In the domain where the arguments of ^ vary, we can for instance 
impose a democratic condition 

\vl\ < ^ (70) 

We remember that p a = —y a , thus condition (70) is a statement about 
individual momenta. From (70) it follows that 

W\ < |ao| + \2y 2 ■ y 3 + -y 2 2 + -y 2 3 \ (71) 

H<N + 2^|^-^| (72) 
Since every y is spacelike, \yA • ]Jb\ < \va\ \vb\- Hence (70) implies 

, ~ ~ !~2 1-2, ^ 

|2y 2 -2/3 + 2^2 + 2^1 < -g" 

|I~2_I~2|<^! 



Therefore 



Now, provided that 



\a-a \<l (73) 
| e — e 1 < ^ (74) 



1 3 



it stems from (73) (74) that a and e remain within the interval (— |, ^). 
In order to realize this situation we are led to restrict the squared-mass 
differences by the condition (75). Then, condition (70) permits to apply 
Proposition I. 

Untill now, we have proposed condition (70) which involves not only the 
relative momenta but also P 2 . Since we consider the positive-energy 
sector of free particles, it is clear that 



^ 2 >E^ = E^c 2 
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For the sake of a simple kinematic interpretation, we have better to 
replace (70) by the stronger condition 

\Pa\ < ^E m 'c 2 (76) 

which is just a little more restrictive and offers the advantage of involv- 
ing only masses and spatial velocities. 

Similarly, in view of (69) it is clear that, in order to fullfill (75), it 
is sufficient to demand 

2 2 

h + ^l <E^ h-d<Ey ( 77 ) 

This approach is well-suited for the equal-mass case and remain useful 
when the mass differences are not too large. 
Example. Two equal masses. 

Assume that m,A = pm\, hence ^^m 2 = (1 + 2p 2 )m\. We find that 

23 25 

(77) is satisfied provided the square- mass ratio satisfies — < p 2 < —. 

It is clear that (76) is a condition on the three-dimensional velocities 
with respect to the rest frame. Although it puts a bound on these 
quantities, it still leaves room for a large class of relativistic motions. 

Example. Three equal masses. In the equal-mass case, m a = m, thus 
both va vanish. We are sure that a, e belong to the safety interval if we 
demand that 

, o, m 2 c 2 . 
\PI\ < — (78) 

Indeed positivity entails that 3m 2 c 2 < P 2 . 

Now what does mean (78) in terms of (Newtonian) velocities ? In 

i ~2 1 2 w2 

the rest frame, for all indices, \p = m ttt^ where w is the 

1 — w z /c z 

cix 1 
Newtonian velocity — . Thus (78) is satisfied provided w 2 /c 2 < -, 

which corresponds to |w| < c/3. Under this limit, say one third of the 
velocity of light, we shall speak of a "moderately relativistic regime". 

For inequal masses, similar results could be derived, but the discussion 
would become a bit complicated. We summarize: 
Proposition II 

In sofar as the mass differences are not too large, we keep the range of 
a, e under control by restrictions on the magnitude of the velocities. If 
in particular we consider three equal masses, velocities under c/3 ensure 
that we can invert our formulas with S(a, e) as in Proposition I. 
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All the quantities involved in condition (70) (resp. (76)) are first inte- 
grals for free particles, thus (70) (resp. (76)) defines an invariant sector 
of the motion. 



4.4 Individuality. New versus old coordinates 

As a result of our transformation of the momenta, it might be puz- 
zling that (beside its dependence on total momentum) each new variable 
q' a depends not only on q a (with the same label a) but also on all q^s 
with b a. This dependence is expressed by the transformation formu- 
las (25). Fortunately, we shall prove that: 
Proposition III 

Beside its dependence on the direction of P, at zeroth order in 1/\P\, 
the variable z' 2 depends only on z 2 (resp. z' 3 depends only on zs). 
Proof 

We develop our formulas in powers of l/\P\ and evaluate z' Aa at lowest 
order. 

dy 

According to (25) we need to compute the coefficients — — . 

dy' 

Let us first prove that 

0(1/\P\), for A^B (79) 



From (42) and (32) it is clear that 

V2 =V2+i P a , Vs = Vt + V P a (80) 

hence 

and similar formulas for dy^/dy'^. We are led to evaluate the derivatives 
of £ (resp. rj). According to (64) it is sufficient to differentiate u and v. 
With help of (50) (51) we get 

P 2 -^r = P a -vt (83) 

(ill 

- ' P «-y%-y% (84) 



dy' 3a 2 
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Let us insert (84) (86) and (83) (85) into the formulas obtained by dif- 
ferentiation of (64). We obtain 



%L-° (1/P) o£r° lllP) <87) 

+ 0(1/ P' z ) (88) 



^ 3P« ln2 . 



dy' 2a 2 P2 

and a similar formula with dr)/dy' 3 . Inserting (87) into (81) we check 



that - . actually is of the order of 1/1 PI, and the same result can be 

derived for dy$/dy' 2 , which altogether proves (79). 

Now we apply formula (25) and take (79) into account. Hence 



A = ^^2 + 0(l/|P|) (89) 

°y 2a 



But in view of (82) (88) we simply have 



o per pa 

1^+0-55-+ 0(1/1^1) 



8y' 2a 2 P2 

So finally 

z' a 2 =z^ + \ ^'p}^ +0(1/\P\) (90) 
and a similar expression in terms of z' 3 ,Z3. In particular we have 

z' 2 = z 2 + 0(l/\P\), z' 3 = z 3 + 0(l/\P\) (91) 



4.5 New form of wave equation 

As seen in Section 4.1, the "difference equations" are (28) or equivalently 
(29) (30). According to (19) the dynamical equation (sum equation) for 
free particles is 

(3j2 m l c2 ~ P 2 ) * = (£> + 6P 2 H) * (92) 
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Of course S must be here considered as a function of y' 2 -P, y^-P, P 2 , VA- 
In view of (56) (60) we can write as well 

where X = S(o~,e) according to (62) (63). We must remember that a, e 
are functions of the new momenta through (49). 

But equation (28) tells that on the mass shell we can replace y' 1A • P 
by vac 2 (thus y' A replaced accordingly, see eqs. (29) (30) ). More- 
over we impose that the total linear momentum has a sharp value k a . 
Let us make this convention that F_ is the expression of any F on the 
momentum-mass shell, namely 

F = subs (y' 1A - P = u A c 2 , P a = k a , F) (94) 

using an obvious notation borrowed from Maple's syntaxis. It is meant 
that y' 1A is as in (22) and we set 

k 2 = M 2 c 2 (95) 

For instance, if we define 

y'a = y'a-(y'-k/k 2 )k a 

we can write y a = y°, therefore 

D = 6[(y' 2 ) 2 + (y' 3 ) 2 + y / 2-y 3 \ 

Moreover (49) yields 



M 2 c 2 u = v 2 c 2 -{y 2 -m + \yl) M 2 c 2 v = ^c 2 -{y 2 -m + \yl) (96) 

It is noteworthy that, in the case of two equal masses e is of the order 
of 1/c 2 , whereas for three equai masses both a and e are 0(1/ c 2 ). 
Taking into account the mass-shell constraints and the sharp value of 
P a we derive the reduced equation 

(3^m 2 -M 2 )c 2 V = (D + 6M 2 c 2 S) V (97) 

Notice that, apart from v 2 , that are fixed parameters, E depends only 
on "ff 2 , y' 3 and M 2 . The only operators involved in (97) are multipli- 
cations by the projections of y' A orthogonal to k, they are essentially 
three-dimensional. Whereas D has a familiar form (just use the rest 
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frame, where Va ■ Vb = ~Ya • Yb) it is not the case for S. Fortunately 
it can be checked that, at least for equal masses, the term c 2 E is in fact 
of the order of 1/c 2 . For this purpose it is convenient to set 



mVs = ^ r (98) 



so we end up with 



(3 ^ m 2 a - M 2 )c 2 = D Y> + ^ (99) 

For three equal masses, T can be expanded in non-negative powers 
of 1/c 2 and it turns out that its zeroth-order piece is biquadratic in ]/ A . 

Proof. It can be easily read off from (49) that in this case u, v thus 
also q_,e are of the order of 1/c 2 . Getting back to system (52) (53) one 
finds that 

| = 2u-v + 0{l/c 4 ) 
r L = 2v-u + 0(l/c i ) 

Inserting into (43) yields 

S = 3(n 2 + v 2 -uv) + 0(l/c 4 ) (100) 

hence 

M 4 c 4 H = T(0) + 0(l/c 2 ) (101) 
T(0) = \ [(y 2 ) 2 + (y 2 ) 2 + 4(y 2 .y 3 ) 2 + 2(y 2 + y 2 ) (fc-fe)-(Sf)(Si)] (102) 

Thus, when all m a = m, the last term in the r.h.s. of (99) can be 
considered as small. 

Free-particle motion is now described only in terms of y 1 and k. 

Imposing by (21) that total linear momentum is diagonal permits, through 
equation (31), to eliminate y' A • k, where the new relative energies cy' A ■ 
k/s/k 2 are conjugate to the new "relative times". 

It is of interest to notice that these new "relative times" are linear 
combinations of the old ones with coefficients that are analytic func- 
tions of the momenta; the reader will check it using (25) (26) (42) and 
(49) (36) (37). 

After reduction, the three-body kinematics has no more degrees of 
freedom than in the non-relativistic problem. But we must keep in mind 
that this picture is valid only in sofar as we can revert to all the initial 
variables, which (at least for equal masses) is ensured for moderately 
relativistic velocities. 

The new variables y' A introduced in this Section will be referred to as 
the reducible variables. 
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5 How to introduce interactions 



We can now consider the system (92) (29) (30) as a starting point for 
introducing mutual interactions. 

To this end, we shall modify the "sum equation" (92) by a term which 
carries interaction, whereas the "difference equations" (29) (30) remain 
untouched. 

Doing so we manage that P remains conserved, and keep assuming that 
its eigenvalue is a timelike vector k; therefore the factorization of ^ given 
by formula (31) remains valid and eliminates two degrees of freedom. 

The interaction potential will be written in closed form in terms of the 
reducible coordinates z' A ,y' B , and all calculations will be carried out 
using these variables. 

Remark: the reducible (momentum) coordinates p' a are re-arranged as 
to form the quantities P and y' A . 

Adding interaction into (92) produces the dynamical equation 

(3 J2 m «c 2 - P 2 )^ = D ^ + (18V + 6P 2 H) * (103) 

Like in the free case, D is given by (17) and S is given by (93) in terms 
of X = S(a,e). 

The "difference equations" remain (29) (30) like previously. Of course, 
V cannot be chosen arbitrarily but it is not difficult to find a general 
admissible form of V such that the dynamical equation (103) is compat- 
ible with (29) (30). Compatibility requires that V commutes with the 
operators in the left-hand sides of (29) (30). For instance the interaction 
potential V may depend on z' 2 , z' 3 and P 2 . 

Naturally V must be Poincare invariant, which is realized by taking a 
function of the various scalar products formed with z' A , y' B , P. 
Demanding that \£ diagonalizes P a with eigenvalue k a , with k ■ k > 0, 
we can in (103) replace y' by y'. 

Taking (29) (30) and (98) into account yields the reduced equation 

6 

(3^m 2 -M 2 ) c 2 y> = £V + 18Z V> + ^2 r ^ (104) 

where the reduced wave function ip depends only on k and on the space 
projections 2/2 > V3 ■ The only operators involved here are the projections 
z' A , y' B . Moreover z' arises in V_ only. 

Comparison with a standard problem of nonrelativistic quantum me- 
chanics becomes more easy in the rest frame, where (2^) 2 = —(z' A ) 2 
and (y' A ) 2 = -(y' A ) 2 , etc. 
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Actually solving (104) differs from a non-relativistic problem by the last 
term, which involves the momenta but does not depend on the shape of 
the interaction (and survives in the free-motion limit). Still this term 
depends on the total squared mass. 
For simplicity, we can consider an interaction such that 

18V = aia U 12 (z' 2 ) + a 23 U 23 (z' 3 -z 2 ) + a 31 U 31 (z' 3 ) (105) 

where a a b are coupling constants and U a b arbitrary (but Poincare in- 
variant) functions. In this model, U± 2 is independent from q' 3 , etc, with 
cyclic permutation; the formal input of our interaction consists in two- 
body potentials. 
So (104) can be written 

(3j2™ 2 a -M 2 )c 2 iP = 

Dip+[a 12 U l2 (^ 2 ) + a 23 Unte-zti + an U 31 {%)] ^ + (106) 
A special case 

U 12 = (z' 2 ) 2 , U 13 = (z' 3 ) 2 , U 31 = (z' 3 - z' 2 ) 2 (107) 
describes a three-boson harmonic oscillator. 

In order to handle equation (104) it is tempting to neglect its last term. 
Invoking the limit of a large total momentum (M 2 —> oo ), as in [12], 
doesnot seem to permit a perturbation treatment. We prefer to consider 
developments in powers of 1/c. 



5.1 Equal masses 

Assuming for simplicity that m a = m, equation (104) becomes 

(9m 2 - M 2 )c 2 ip = Ihp + 18V tp + T i/j (108) 

It will be considered as an eigenvalue problem for A by setting 6A = 
[M 2 — 9m 2 )c 2 . As all masses are equal, thus T = T(0) + 0(l/c 2 ). At 
first order in 1/c 2 we can replace T by T(0) and M 2 by 9m 2 in (108). 
Using the rest frame we obtain 

= [yl + y§ + y 2 ■ y 3 - 3Z] V - ^2 r (°) ^ (109) 

Neglecting the last term yields the nonrelativistic limit (divide by m 
and remember that in our formulas, V has dimension of P 2 ). 
Taking into acount the contribution of T(0) permits to calculate the 
first relativistic correction. 
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5.2 Cluster behaviour 



As pointed out by Sazdjian [10], in any formulation of the dynamics 
which makes explicit reference to total momentum, it is difficult to dis- 
cuss cluster separability. But it is reasonable to demand that the reduced 
equation be in a sense separable, in order to ensure a factorization of 
the internal wave function when there are noninteracting clusters. 
With this requirement in mind, we can already observe that the poten- 
tial (105) is formally separable in terms of the variables z'. 
But the interpretation of each U a b as a two-body term runs into a com- 
plication: there is no evidence that the variable z' A , exactly matches 
the cluster of particles {1A}. A similar remark arises concerning the 
matching of z' 3 — z' 2 with cluster {23}. 

The physical interpretation of the new configuration variables z' A is not 
straightforward; they are relative variables since they commute with P, 
but they suffer from this complication that the transformation formulas 
(25) mix Z2 with zj,. Similarly (beside its dependence on total momen- 
tum) each new variable q' a depends not only on q a (with the same label 
a) but also on all q^s with b 7^ a. 

But we can consider (91) on the momentum-mass shell. At least for 
three equal masses, the only occurence of the velocity of light is through 
the product Mc, so we obtain from (91) 

% = z 2 + 0(l/\Mc\), zf 3 = z 3 + 0{l/\Mc\) (110) 

and, of course, z^ — z^ = z% — Z2 + 0(l/\Mc\). Thus, at leading order, the 
variables z! A and za still coincide; so the potentials U a b in the reduced 
equation (106) can be approximately considered as two-body terms. 



6 Concluding remarks 

As a first step, we succeeded in constructing three mass-shell constraints 
describing the free motion of three scalar particles. In contrast to the 
KG system, these new wave equations permit to eliminate two degrees of 
freedom and get reduced to a covariant equation with three-dimensional 
arguments. 

Our approach rests on a transformation of the momenta involved in 
the original KG system. In contradistinction to Sazdjian's proposal and 
the homographic relations that approximate it (eq (13) of ref. [10], eq 
(4.15) of ref. [11]), our transformation from the old momenta to the 
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new ones is explicitly given by simple quadratic formulas. Inversion of 
these formulas is a fourth degree algebraic problem which could be (in 
principle) discussed and solved in closed form; due to its complexity, 
approximate developments are more efficient in practical calculations. 
We used a couple of identities that are specific of the three-body case; 
thus an extension of the present work to n > 3 is by no means straight- 
forward! 

In the present state of the art, equivalence of the new equation (97) 
with the sum of the original KG equations is ensured at least in a large 
sector characterized by positive energies and conditions that involve the 
masses of the particles. When the masses are not too different one from 
another (and in particular for equal masses), these conditions amount 
to impose a bound on the velocities; but this bound is still high enough 
to allow for the description of a relativistic regime. 
The case of very large velocities requires further investigations. We gave 
here sufficient conditions for an invertible transformation; it remains 
possible that a more detailed discussion enlarges the present results. 
This analysis of free-body kinematics provides us with a solid ground. 

In a second step, we introduced interaction in the "sum equation". The 
model obtained by this procedure respects Poincare invariance. It re- 
mains covariantly reducible to a wave equation with three-dimensional 
arguments; free motion is recovered in the absence of interaction term. 

The interaction term V is formally cluster separable; actually for- 
mula (105) is an ansatz which permits to combine two-body interactions 
without spoiling the compatibility of the mass-shell constraints. True 
separability (in terms of the original individual particle coordinates) is 
recovered only in the large-total-mass limit. 

The two-body input of our model can be either phenomenological or 
motivated by consideration of field theory. 

When the three masses are equal, the velocity of light arises in S 
through the product M 2 c 2 , which facilitates the expansion in powers of 
1/c 2 . At the first order, the reduced equation is similar to a familiar 
Schroedinger equation supplemented with a perturbation; insofar as the 
interaction is not explicitly energy dependent (or if this dependence is 
of higher order) one is left with a conventional eigenvalue problem. 
This situation provides a basis for eventually undertaking the study of 
cases where the mass differences are not zero but still remain relatively 
small. 

In the hope of applications to three-quark or three-nucleons systems, we 
plan an extension of the formalism to particles with spin. The contact 
with more elaborated (but more complicated) theories, such as QED 
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and QCD, will be discussed in a future work. 
Appendix 1 

In non-relativistic classical mechanics, canonical transformations are 
symplectic diffeomorphisms of phase space. In general they do mix 
the q's and the p's. 

But a point transformation (in configuration space) simply transforms 
the q's among themselves, say q' = f(q). Then invariance of the sym- 
plectic form [13] fully determines the p"s in terms of the q 's and the p's. 
When configuration space is flat, the q and p variables play symmetric 
roles in the general formulas of analytic mechanics, so there is no diffi- 
culty in defining as well point transformations in momentum space (but 
this possibility is not usually considered in textbooks). In this case, one 
transforms the momenta among themselves, and one further determines 
the new variables q' in terms of g's and p's through the requirement that 
the complete transformation law is canonical. 

In the position (resp. momentum) representation of quantum mechan- 
ics, a quantum analog of point transformations in configuration (resp. 
momentum) space can be generated by an invertible transformation of 
the arguments of the wave function. This transformation among c- 
numbers obviously induces a transformation among the multiplicative 
operators they define. 

Appendix 2. 

The polynomial R <T {X) has an obvious double root X = 2/3 inde- 
pendent of a, and provided a > — 1, two other real roots given by (59) 
but, as noticed above, the root X~ falls outside the admissible interval. 
All the curves Z = R a (X) are tangent to the X axis at a fixed point 
X = 2/3. For a > — 1 and e small enough, the curve representing 
R(X) is four times cut by the straight line Z = e 2 . In the limit when 
e vanishes, two points of this intersection form the contact with the X 
axis, and the other ones respectively reduce to X + and X~. 
For a = ^ we find that X + ~ 0.15, which is admissible in the sense of 
(45) (46), and Ri (X) has a local maximum at X = 1/3. This maximum 

is exactly \. For a < ^, the local maximum exceeds \. Making a to 
decrease we obtain lower values of X + (which vanishes with a). 
For a = — |, we obtain exactly X + = — |, and going down further is 
excluded in view of (46). 

Taking a in the open interval (— |, |) and X restricted by — | < X < |, 
it turns out that, provided e 2 < |, each curve Z = R a (X) has two 
points in common with the straight line Z = e 2 (other possible points 
correspond to X outside the interval we consider). For vanishing e, 
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one of them has its horizontal coordinate going to coincide with X + , 
while the other point goes to the fixed contact point X = 2/3, Z = 0. 
This analysis shows that, with our restrictions, the 4th dgree equation 
R a (X) = e 2 has two real solution, but only one of them reduces to X + 
in the limit where e vanishes. 
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